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Inhomogeneous metal enrichment at z ~ 1.9: the Lyman limit 
systems in the spectrum of the HDF-S quasar* 
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Abstract. We present a detailed analysis of three metal absorption systems observed in the spectrum of the HDF- 
South quasar J2233-606 (z cm = 2.238), taking advantage of new VLT-UVES high resolution data (R = 45000, 
S/N= 40 - 60, AA3050 - 10000 A). Three main components, spanning about 300 km s , can be individuated 
in the Lyman limit system at z ~ 1.92. They show a surprisingly large variation in metallicities, respectively 
~ 1/500, 1/8 and 1/100 solar. The large value found for the second component at z ~ 1.9259, suggests that the 
line of sight crosses a star-forming region. In addition, there is a definite correlation between velocity position and 
ionisation state in this component, which we interpret as a possible signature of an expanding Hn region. The 
systems at z ~ 1.94 and z ~ 1.87 have also high metallicity, ~ 1/4 and 1/3 solar. We find that photoionisation 
and collisional ionisation are equal alternatives to explain the high excitation phase revealed by O VI absorption, 
seen in these two systems. From the width of the Siiv, Civ, Sim and Cm lines in the system at z ~ 1.87, we 
can estimate the temperature of the gas to be logT < 4.7, excluding collisional ionisation. Finally, we compute 
the Siiv/Civ ratio for all Voigt profile components in a sample of logiV(Civ) > 14 systems at z < 2. The values 
show a dispersion of more than an order of magnitude and most of them are much larger than what is observed for 
weaker systems. This is probably an indication that high column density systems preferably originate in galactic 
halos and are mostly influenced by local ionising sources. 
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1. Introduction 
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High Hi column density systems (JV(Hi) > 10 
and especially Lyman limit systems (AT (Hi) > 10 
cm~ 2 , hereafter LLS) arise in dense environments such as 
halos of large galaxies or the densest regions of filaments ionisation state and abundances (e.g. [Petitjean ct al. 1992 



Although the ionisation corrections are large in LLS, 
the large number of associated metal lines can be used 
to constrain ionisation models and to estimate metallicity, 



linking the galaxies (e.g. Katz et al. 1996; Gardner et al. [Bergeron et al. 1994|; [Kirkman k Tytlcr 1999j; |Kohler et al 



1997] ) ■ |Indeed, at redshifts z < 1, galaxies associated with 1999| ; |Prochaska k Buries 1999|; |Chen k Prochaska 2000|) 



LLS (detected by Mg n absorption with equivalent width 
w T > 0.3 A) are routinely identified revealing the presence 
of gaseous halos with radius larger than w 40 h^ 1 kpc 
([Bergeron k Boissc 199l|; |Steidel et al. 1994|; phurchill 



UVES, the UV and Visible Echelle Spectrograph, mounted on 
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In this paper, we present the study of three metal 
systems at redshifts z a b s ~ 1-87, 1.92, and 1.94 seen in 
the spectrum of the HDF-South quasar J2233-606 (z cm = 
2.238, B ~ 17.5; see |Savaglio 199S| ; |5avaglio et al. 1999Q . 
The systems are associated with a Lyman limit break 
observed at ~ 2683 A ( [gealey et al. 19"98 ). Prochaska 



Bas ed on material eollcetcd during commioaioning time o 



the Kucvcn ESQ telescope operated on Cerro Paranal (Chile' 



and^vitli^lM^ASA^ESA^M 



at the Space Telescope Science Institute, which is operated by 



the Association of Universities of Research in Astronomy, Inc 



under NASA contract NAS 5-26555 



k Buries (1999] , PB99) investigated the two systems at 
z ~ 1.92 and z ~ 1.94 in detail. They derived logAT(Hi) 
= 17.15 ± 0.02 and 16.33 ± 0.04 at z ~ 1.92 and 1.94 re- 
spectively. They found the systems have solar abundance 
pattern, but significantly different metallicites: ~ 1/50 so- 
lar and less than w 1 /200 solar for the two components at 
z - 1.92 and w 40 % solar at z ~ 1.94. 
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Table 1. Ionic column densities for the z ~ 1.92 system. 
An (f) indicates that in the fit the parameter has been 
fixed to the value obtained for the previous absorption 
line 



Comp. 


z 


b 

(km s- 1 ) 


Ion 


log AT 
(cm" 2 ) 


1. . . 


1.92543" 


11 ± 5 


Si IV 


11.7 ±0.2 








Sim 


11.2 ±0.2 








Civ 


12.2 ±0.2 


2. . . 


1.925719 6 


14 ±1 


Si iv 


12.23 ±0.05 








Sim 


12.61 ±0.04 








Civ 


13.07 ±0.03 




(f) 


(f) 


Nv 


< 12.75 ±0.08 




(f) 


(f) 


Nm- 


13.5 ±0.2 




1.925718 6 


18 ±11 


di 


12.4 ±0.2 


3... 


1.925950 c 


7.6 ±0.6 


Sin 


12.26 ±0.03 








Cn 


12.91 ±0.05 








Mgn 


11.72 ±0.08 




1.925989 c 


6.7 ±0.2 


Si iv 


12.72 ±0.01 








Sim 


12.76 ±0.09 








Civ 


13.41 ±0.01 




(f) 


(f) 


Nv 


< 12.4 ±0.1 




(f) 


(f) 


Nm 


13.6 ±0.2 


4. . . 


1.927747 c 


7.4 ±0.8 


Sim 


12.32 ±0.03 


5... 


1.927929 c 


5.5 ±0.8 


Sim 


12.14 ±0.04 
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a Error 30 x 10" 



Error 7 x 10" 
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Fig. 1. Main ionic transitions observed for the system at 
z ~ 1.92. The dotted lines mark the position of the clouds 
discus aod in the text. From loft to right: cloud 1 at v — 43- 
km {z — 1.0255) with motallicity ° 1 1/500 polar; cloud 
2 at v = km s" 1 (z = 1.9259) with metallicity - 1/8 
solar and cloud 3 at v ~ ±184 km s -1 (z = 1.92577) with 
metallicity ~ 1/100 solar 

Here we take advantage of new complementary data 
of high resolution and high signal-to-noise ratio (S/N) to 
further constrain metallicity inhomogeneities in these sys- 
tems. The paper is structured as follows: in Sect. 2 we 
report observational details and the line fitting procedure 
for the systems; in Sect. 3 we investigate the possible ionis- 
ing mechanisms and the determination of metallicity and 
abundances. In Sect. 4, the observed behaviour of the ratio 
Siiv/Civ is described; finally, we report our conclusions 
in Sect. 5. 

2. Comments on individual systems 

New data on the HDF-South quasar have been obtained in 
October 1999 during the commissioning of UVES, the UV 
and Visual Echelle Spectrograph, mounted on the VLT 
Kucyen ESO telescope at Paranal (Chile). The spectrum 



is of high resolution (R = 45000) and high signal-to-noise 
ratio (S/N = 40 — 60 per resolution element) and covers 
the spectral range AA3050 — 10000 A. The data reduc- 



tion and spectrum analysis are reported in Cristiani & 



D'Odorico (2000| ). We use these data together with the 
echelle spectrum (AA2275 — 3118 A) of high resolution 
(R = 30000) obtained with the STIS instrument on board 
HST qSavaglio 1998| ). 

When not otherwise stated, low-ionisation, weak lines 
of different species have been fitted together with the same 
number of sub-components, each with the same redshift 
and Doppler width. The same is true for high ionisation, 
complex and/or saturated lines. Absorption lines detected 
in the STIS spectrum (mainly N III, C ill) are fitted using 
the same components used to fit Si IV. Voigt profile fit- 
ting is obtained in the context Lyman of MIDAS the ESO 
reduction package (Fontana & Ballester 199!:). 



2.1. The system at z ~ 1.92 

The strong, well defined component at z ~ 1.92599 (see 
Fig. |l| and Table |l|) in the high-ionisation line profiles 
(Civ, Sim, Siiv) is displaced by 4 ± 0.4 km s^ 1 with re- 
spect to the corresponding component at z ~ 1.92595 seen 
in the low ionisation lines (Cli, Sin) as shown in Fig. [l| 
The number of lines available together with the very good 
quality of the data give confidence that this shift is real. 
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Fig. 2. Main ionic transitions observed for the system at 
z ~ 1.94. The origin of the velocity axis is fixed at z — 
1.942616 and the dotted lines indicate the positions of 
components 7 and 8 in Tabic ^ 

The simplicity of the velocity profiles suggests that this re- 
flects the internal structure of an H n region flow in which 
kinematics and ionisation state are well correlated (e.g. 
Hcnney fc O'Dell 1999| ; |R.auch et al. 1999| ). 

The simultaneous fit of the Hi Lya, Ly/3, Lye and Ly- 
8 lines results in two main components at z ~ 1.92553 
and z ~ 1.92775 with logAT(Hl)= 17.1 ± 0.2 and 16.7 ± 
0.3 respectively. Thus the two stronger Hi absorptions 
are shifted relative to the main C n and Si n component 
observed at z = 1.92595 by ~ -43 km s _1 and +184 km 
s _1 respectively. The mismatch in the velocity positions of 
the strongest H I and C n absorptions is apparent in Fig. |l|, 
where the origin of the velocity axis is fixed at the position 
of the C II absorption. The upper limit to the H I column 
density at the redshift of the main metal component is 
logTV(Hi) < 16. This is derived by adopting the Doppler 
parameter as measured from the C n line. It can be seen 
as well that at velocity position +184 km s _1 , the only 
metal transition detected is Si III. The 3 a upper limit on 
the corresponding Civ absorption is logiV(Civ) < 11.9. 



Fig. 3. Main ionic transitions for the system at z ~ 1.87. 
The origin of the velocity axis is fixed at z = 1.87008 
corresponding to component 12 analysed in the text 




Relative Velocity (k: 



Fig. 4. Superposition of the absorption lines due to Si n 
A 1260, Cn A 1334 (solid line) at z - 1.92595 and Siiv 
A 1393, Civ A 1548 (short-dashed line) at z - 1.92599. 
The dot-dashed lines indicate the center of the two groups 
of lines separated by 4 km s _1 
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Table 2. Ionic column densities for the z ~ 1.94 system 







b 




log iV 


Comp. 


z 


(km s" 1 ) 


Ion 


(cm" 2 ) 


1. . . 


1. 94031° 


12 ±1 


Civ 


12.67 ±0.04 




1.940258 6 


19 ±1 


Ovi 


13.9 ±0.1 








Hi 


13.28 ±0.03 


2. . . 


1.940730 6 


7.0 ±0.4 


Civ 


12.97 ±0.02 






22 ±3 


Ovi 


13.9 ±0.2 








Hi 


13.53 ±0.07 


3... 


1.941006 6 


10 ±1 


Civ 


12.78 ± 0.03 




1.94107 c 


22 ±4 


Ovi 


14.1 ±0.1 








Hi 


13.2 ±0.1 


4... 


1.941791 6 


5.8 ±0.6 


Civ 


12.98 ±0.05 


5. . . 


1.941975 6 


14.3 ±0.7 


Si iv 


12.78 ±0.03 








Sim 


12.70 ± 0.03 




(f) 


16.4 ±0.6 


Civ 


13.89 ±0.01 




1.94182 c 


20 ±4 


Cii 


12.9 ±0.1 




(f) 


91 ± 25 


Hi 


13.80 ±0.07 


6. . . 


1.942001 6 


3.7 ±0.5 


Si iv 


12.66 ±0.04 








Sim 


12.4 ±0.1 




(f) 


5±1 


Civ 


12.84 ±0.08 




1.942011 6 


9±2 


Cii 


12.8 ±0.1 








Sin 


11.79 ±0.08 




(f) 


23.9 ±0.5 


Hi 


15.60 ±0.06 


7. . . 


1.942429 6 


10.3 ± 0.4 


Si iv 


13.34 ±0.02 








Sim 


13.20 ±0.04 




(f) 


(f) 


Nm 


13.7 ±0.1 




(f) 


10.1 ± 0.4 


Civ 


14.01 ±0.03 




1.942484 6 


17.1 ±0.6 


Cii 


13.62 ±0.02 








Sin 


12.59 ±0.03 








Mgn 


12.54 ±0.04 








Aim 


11.92 ±0.05 








Aln 


11.39 ±0.05 




(f) 


26 ±2 


Hi 


16.1 ±0.1 


8... 


1.942618 6 


9.1 ±0.03 


Si iv 


13.39 ±0.02 








Sim 


13.46 ± 0.05 




(f) 


(f) 


Nm 


13.8 ±0.1 




(f) 


14 ± 1 


Civ 


14.1 ±0.1 




1.942614 6 


5.2 ±0.3 


Cii 


13.52 ±0.03 








Sin 


12.53 ±0.03 








Mgn 


12.37 ±0.05 








Aim 


11.79 ±0.05 








Aln 


11.21 ±0.06 




(f) 


15 ±4 


Hi 


15.8 ±0.2 


9... 


1.94269 c 


18 ±1 


Civ 


13.8 ±0.2 


Error 11 


-20x 10" 6 , b 


Error < 10x10" 6 , c 


Error > 30 x 10 



2.2. The system at z ~ 1.94 

Absorption profiles for this system are shown in Fig. ^|. 

Si IV and Si in absorptions are fitted with the same 
components and Doppler parameters. Civ is saturated; 
we therefore fix the components at the redshifts deter- 
mined for Si iv and then add weak components to reach a 
good fit. Sin A 1193, A 1260, A 1304, A 1526, Cn A 1334, 
Mgn A 2803, Aln A 1670 and Aim A 1854 are fitted to- 




Relative Velocity (km s ') 



Fig. 5. O vi AA 1031, 1037 for the system at z ~ 1.94. The 
velocity axis is centered at z = 1.941975. The dotted lines 
show the position of the main components in C IV 



gether (same redshift and same Doppler parameter) . Mg n 
A 2796 is blended with a telluric absorption and is not 
considered in the fit. Results are given in Table |. Most 
components shows significant differences in velocity (up to 
~ 25 km s _1 ) between high-ionisation and low-ionisation 
transitions, but blending is likely the main cause of this 
discrepancies. 

Associated O VI is detected in the STIS portion of the 
spectrum (see Fig. ||). Although the S/N ratio is not ex- 
cellent, the Ovi and Civ profiles are correlated with in 
particular three components (1, 2, 3 in Table |2|) well de- 
tected in both species at z - 1.94026, 1.94073, 1.94107 
(v = —175, —127 and —92 km s -1 respectively in Fig. |^). 

The Hi column densities quoted in Table || are ob- 
tained by fitting together Hi Lya, Ly7, Lye, Ly-10, Ly-12 
at the redshifts fixed by the Cii components. Nv is not 
observed; the 3 a upper limit on the column density is 
logN(Nv) < 12.15. 

2.3. The system at z ~ 1.87 

This system has a complex structure spanning around 600 
km s _1 , with 11 (16) components in Si IV (Civ). Mgn is 
observed only at the redshift of the two central compo- 
nents (11 and 12 in Table Q). Unfortunately, important 
transitions like Cii A 1334, Sill A 1260, Sim A 1206 and 
Cm A 977 are badly affected by blending (see Fig. ||). 

The Si iv doublet is fitted by itself and the redshifts 
of the components are used to fit Si III. C iv is fitted by 
itself because more components with respect to Si iv are 
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Fig. 6. O vi AA 1031, 1037 for the system at z ~ 1.87. The 
velocity axis is centered at z = 1.8695. The dotted lines 
show the position of the main components in C iv 

necessary to reach an acceptable fit. Table ^ displays the 
fitting parameters for the 17 components. 

Associated O VI is detected in the STIS portion of the 
spectrum (see Fig. g), low S/N ratio and partial blending 
with intervening H I lines prevent to carry on a best fit 
procedure in MIDAS. An indicative upper limit to the 
column densities in each component is estimated t o be 
13.7 using t he interactive fitting program Xvoigt ( Mar 



& Bail 3^ 1995| ). It is apparent from Fig. gjthat, although 



the C iv and O vi profiles are correlated, the distribution 
of the O VI phase is smoother than that of the C iv phase. 
Nv is not detected: the 3cr limit on the column density 
is: logiV(Nv)< 12.3 for z < 1.87 and logiV(Nv)< 12.6 
for z > 1.87. 

Due to the heavy saturation of Lycn and Ly/3 absorp- 
tions and to the complexity of the system, it is not possi- 
ble to reliably constrain Hi column densities for all the 
components. For components 2, 4 and 16 we estimate 
logA(Hi)~ 15.2 ± 0.2, ~ 14.9 ± 0.1 and ~ 15.4 ± 0.1 
respectively from the simultaneous fit of Hi Lya, Ly/3, 
Ly7, and Ly-6. 

3. lonisation models 

The nature of the ionisation mechanisms is investigated 
comparing the observed column densities with those pre- 
dicted by photoionisation models. For this we use the 
Cloudy software package ( Fcrland 1997 ) assuming ioni- 
sation equilibrium. In the models we adopt: (1) a total 
hydrogen density n(H)= 10~ 25 cm -3 , (2) a plane paral- 
lel geometry for the gas cloud, (3) the solar abundance 



pattern, and two ionising spectra (4a) the Haardt-Madau 
( [Haardt fc Madau 1996] ; [Madau et al. 1999| , hereafter HM) 
EUVB spectrum with galaxy contribution or (4b) a com- 
posite (hereafter SSP) of the spectrum produced by a sim- 
ple stellar population with 20 % solar metallicity, age 0.1 
Gyr (Chariot private communication) plus a power law, 
f(y) oc v , which increases the flux at energies greater 
than the Hen ionisation break (54.4 eV) by a factor of 
~ 3 (|5chaerer fc de Koter 19971), 



U 



The ionisation parameter, U, defined by 

_ 0912 _ 47TJ 912 _ 2 x 10 _ 5 J912/10" 215 



cn H 



(1) 



hcna Tin/ (lcm 3 ) 

(where J912 is the intensity of the EUVB at 1 ryd) and 
the metallicity of the system are derived by matching the 
model output column densities with the observed values. 

3.1. The system at z ~ 1.92 

In Sect. 2.1 we have singled out three clouds at z ~ 1.9255, 
1.9259, 1.9277 with log AT (Hi) ~ 17.1, <16, 16.7 respec- 
tively, corresponding to the v = —43 km s _1 component 
which is the main component in Hi (cloud 1), the v = 
km s _1 component which is the main component in C11 
and Sin (cloud 2) and the +184 km s _1 component which 
is detected in Hi and Si in only (cloud 3). At variance with 
PB99, we anlyse separately cloud 1 and cloud 2. As there 
is an uncertainty of ~ 20 km s _1 on the redshift determi- 
nation of cloud 1, we use for it the metal column densities 
of the slightly shifted component 2 at z ~ 1.92572 (see 
Tab. 0). Thus, the metallicity we determine for this cloud 
could be an upper limit. Finally, the model for cloud 3 
takes into account the possible Si ill absorption line (com- 
ponent 4 in Tab. |l|) which was not observed by PB99. 

Because of the large number of parameters (ionisa- 
tion parameter, shape of the ionising flux, metallicities), 
the metallicity determination is most often degenerate. 
However, when the number of transitions is large and the 
column densities well determined, the mean metallicity 
and ionisation parameter can be estimated with reason- 
able uncertainties. 

The observed column densities are reproduced by 
a SSP spectrum and a slight overabundance of silicon 
([Si/H] ~ 0.15). The higher quality data allow us to re- 
ject the HM spectrum model which was considered vi- 
able in the analysis by PB99. Indeed, the observed column 
density pattern (Cn, Civ, Sin, Sim, Siiv) of cloud 2 is 
not reproduced by the HM spectrum for any combination 
of metallicity and ionisation parameter, even allowing for 
an over-solar abundance of silicon. Fig. [f] shows the pre- 
dicted column densities as a function of U for cloud 2 mod- 
els. The column densities obtained in our best model are 
given in Table [| The computed values for silicon ions can 
still be reconciled with the observed profiles. The discrep- 
ancy between observed and predicted N V and N ill column 
densities, discussed in PB99, is no longer present in our 
analysis. It is apparent from Fig. that the Al 11 column 
density is predicted larger than the upper limit derived 
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Table 3. Observed and predicted ionic column densities 
for cloud 1 in the system at z ~ 1.92 (see text) and com- 
ponent 8 in the system at z ~ 1.94 





cloud 1 (z ~ 


1.92595) 


comp. 8 (z ~ 


1.942618) 


Ion 


log iV b s 


log N pTed 


log 7V obs 


log iVp re d 


ClI 


12.91 ±0.05 


12.94 


13.50 ±0.03 


13.5 


Civ 


13.41 ±0.01 


13.41 


14.1 ±0.1 


13.8 


Niii 


13.6 ±0.2 


13.48 


13.8 ±0.1 


14.2 


Nv 


< 12.4 ±0.1 


11.55 


< 12.15 


11.56 


Mgn 


11.72 ±0.08 


11.67 


12.37 ±0.05 


12.39 


Aim 


< 11.19 


10.92 


11.79 ±0.05 


12.09 


Aln 


< 10.66 


11.25 


11.21 ±0.06 


11.82 


Sin 


12.26 ±0.03 


12.16 


12.59 ±0.04 


12.59 


Sim 


12.8 ±0.1 


13.1 


13.46 ±0.05 


13.81 


Si iv 


12.72 ±0.01 


12.64 


13.39 ±0.02 


13.39 



Fig. 7. Ionic column densities computed by Cloudy 
(Ferland, 1997) for \ogN m = 16, [X/H] = -0.87, [Si/H] = 
0.15 and an SSP ionising spectrum with 20 % solar metal- 
licity (see text) vs. a modified ionisation parameter, U, as 
defined in the text. The symbols on the right side of the 
plot and the dashed orizontal lines indicate the observed 
column density values for cloud 2. The symbols sifhted to 
the right are upper limit measurements. The vertical dot- 
ted line indicates the best value of U determined from the 
observed Civ, ClI, Siiv, Siniand Sin column densities 

from observations. This behaviour is even more apparent 
in Fig. U both for the column densities of Al in and Al II 
and for their ratio which is predicted ~ 0.3 dex lower 
than observed. The underabundance of aluminium ob- 
served in galactic low metallicity stars could reconcile the 
observed and computed results. On the other hand, mag- 
nesium does not seem to be overabundant as would be ex- 
pected from the same halo star observations. Furthermore, 
the recombination coefficients used to compute the alu- 
minium ionisation e quilibrium are probab ly questionable 
as already noted by Petitjean et al. (1994 ). 

For the other clouds, we derive that cloud 1 has similar 
ionisation parameter as cloud 2, \ogU ~ —2.7, —2.8. The 
metallicities of the two clouds differ by two order of magni- 
tudes: [X/H] ~ —2.7 and —0.9 respectively. The computa- 
tion for cloud 3 is based solely on the Si in column density 
and on the upper limits on the main transitions of silicon 
and carbon. The gas is probably in a low ionisation state, 
—3.3 < log!/ < — 3.1, and has a low metallicity —2.3 < 
[X/H]~< -2.0. 

From this we derive that most probably the three 
clouds are ionised by a local source of ionisation. This 
conclusion is strengthened by the observed correlation be- 



tween kinematics and ionisation state in the component at 
z ~ 1.9259 (cloud 2). This together with the observations 
of variations in the metal content by about two orders of 
magnitude over a velocity interval of ~ 50 km s _1 , sug- 
gests that the line of sight crosses a region of intense star- 
formation activity. It would be most interesting to probe 
this by searching the field for emission line objects at this 
redshift. 



3.2. The system at z ~ 1.94 

Photoionisation models have been constructed for the 
two strong components 7 and 8 (see Table ^) both with 
logA^(Hi) « 16. We have used as constraints the column 
densities of Sin, Cn, and Siiv. We do not use Civ and 
Si III because the lines are saturated. 

We cannot find any satisfactory fit to the data for a 
SSP ionising spectrum, except allowing for a Si/C ratio 
smaller than solar which is probably unreasonable. On the 
other hand, the HM spectrum accounts for the observed 
column densities in both components for log U ~ —2.3 and 
—0.7 < [X/H] < —0.5. The uncertainty on the metallicity 
is mostly due to the uncertainty on the H I column density 
determination. The model requires also an overabundance 
of silicon, [Si/H] ~ 0.2. In Fig. || the computed column 
densities for all the observed ions as a function of U are 
reported in the case of component 8 and the values for 
the best fit parameters are in Table ||. It can be seen that 
most of the column densities are explained within a factor 
of two which is probably acceptable given the uncertainties 
on the column densities. 

As previously noted, there is a discrepancy between the 
observed and computed Alii and Aim column densities. 
Niii is slightly overpredicted. Nitrogen could be under- 
abundant which would not be surprising in the hypothe- 
sis of secondary origin of this element. A firm conclusion 
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Fig. 8. Results from a series of Cloudy (Ferland, 1997) 
calculations assuming an overabundance [Si/H] ~ 0.2, 
a metallicity [X/H]=-0.5, log7V ffl = 15.8 and an HM 
ionising spectrum including the contribution of galaxies. 
Plotted are the calculated ionic column densities vs. a 
modified ionisation parameter, U, as defined in the text. 
The symbols on the right side of the plot and the dashed 
orizontal lines indicate the observed column density values 
for the component of the metal system at z = 1.942614, 
symbols shifted to the right give upper limit measure- 
ments. The vertical dotted line indicates the best value 
of U determined from the observed Cn, Siiv and Sin 
column densities 



cannot be reached as the N in column density is not well 
determined. 

The three weak components 1, 2, and 3 show ab- 
sorption due to Ovi, Civ and Hi. The observed ratios 
log (Civ/Ovi), log (Civ/Hi) are similar for the three 
components, average values are ~ —1.1 ±0.2 and —0.5 ±0.1 
respectively. We investigate both the possibilities of photo- 
or collisional ionisation. 

The two adopted ionising spectra can reproduce 
the observed column density pattern in the context 
of photoionisation. But Nv is always overpredicted: 
logiV(Nv)~ 13 instead of <12.15 (3<r) as observed. In 
the case of the HM spectrum the parameters vary in the 
ranges -0.15 < logU < 0.2 and -1.4 < [X/H] < - 1.1 
for the three components. While for the SSP spectrum 
-0.86 < log U< - 0.55 and -0.52 < [X/H] < -0.22. 

In the hypothesis of collisional ionisation, the observed 
column density ratio log (C iv/Ovi) is consistent with a 
temperature logT = 5.35 ( [Sutherland fc Dopita 1993| ). 
The prediction for the other ratios are: log (Civ/Hi)~ 



—0.03, log (O vi/Hi)~ 1.30 for solar metallicities. In order 
to match the observed ratios, the metallicity should be de- 
creased to ~ 1/3 solar. Thus, if we assume collisional ion- 
isation for these components and a solar abundance pat- 
tern, we get a metallicity comparable to the one predicted 
by Cloudy with the SSP spectrum. At this temperature, 
collisional ionisation predicts also: log (Nv/Civ)~ 0.6 
while —0.5 is observed. Therefore in the framework of 
these models, Nv observations can be explained only if 
nitrogen is deficient by an order of magnitude compared 
to carbon. 

In conclusion, the gas in this system is likely of quite 
high metallicity (larger than 0.1 solar) for this redshift. 
The ionisation state can be explained by photoionisation 
by a HM type spectrum. In this case however, metallicity 
is three times less in the high-ionisation phase (compo- 
nents 1, 2 and 3) compared to the low- ionisation phase 
(components 7 and 8). The model can be reconciled with 
similar metallicities in both phases in case the high ioni- 
sation phase is predominantly ionised by a local source of 
photons or if it is collisionally ionised. 

Both the photo- and collisional ionisation models pre- 
dict that nitrogen is underabundant with respect to solar 
[N/C] < - 0.5. 



3.3. The system at z ~ 1.87 

We first determine an upper limit on the temperature of 
the gas from the measured Doppler parameter, in the hy- 
pothesis of pure thermal broadening. We choose three nar- 
row, resolved components (2, 4 and 17; see Table || and 
Fig. ||) with detected Siiv, Sim, Civ and, the last one, 
C in absorptions; the observed Doppler parameter are con- 
sistent with a temperature logT < 4.7. In the hypothe- 
sis of collisional ionisation, at this temperature most of 
the carbon would be in Cn and Cm ions, Ovi would 
not be present, and silicon would be mainly in Sim and 
Siiv. From the ionic fractions calculated by [Sutherland 
|fc Dopita (1993| ), at logT = 4.6 log (Sim/Siiv) ~ 0.8, 
and at logT = 4.8 log (Sim/Siiv) ~ 0.16 and log 
(Cm/Civ)~ 1.74. This excludes the possibility that the 
gas phase giving rise to the observed silicon and carbon 
transitions is collisionally ionised. 

However, Ovi could still arise from a collisionally 
ionised phase. Indeed, the co-existence of these two ionis- 
ing mechanisms has been explained by e.g. the presence of 
a tail of shock-heated gas at high temperature away from 
the curve of equilibrium ( Hachnelt ct al. 1996] ) or a two 
phase halo where clouds photoionised by the EUVB at an 
equilibrium temperature of ~ 10 4 K and giving origin to 
the relatively narrow lines of Civ, Siiv, etc., are in pres- 
sure e quilibrium with the hotter halo gas in which O v i 
arises flPctitjean et al. 1992| ; [Mo fc Miralda-Escude 1996|) . 
The latter scenario predicts broad, shallow absorption for 
O vi and C iv originating in the hot gas. The low S/N 
ratio STIS spectrum marginally shows that these lines are 
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log 10 (Mgll/CIV) log 10 (CII/CIV) 




log 1D (Silll/SilV) 

Fig. 10. Results from a series of Cloudy calculations 
(Ferland, 1997). Curves represent Sin/Siiv vs. Sim/Siiv 
column density ratio as the ionisation parameter varies in 
the range —4 < \ogU < — 1. The different curves describe 
the same models as in Fig. | for an HM ionising spectrum 
(panel a) and an SSP spectrum (panel b) . The superposed 
point is the corresponding observed ratio 



Fig. 9. Results from a series of Cloudy calculations 
(Ferland, 1997). Curves represents Siiv/Civ column den- 
sity ratio vs. Mg n/C iv (panels a-d) and C Ii/C IV (panels 
e-h) ratio as the ionisation parameter varies in the range 
—4 < \ogU < — 1. Different curves are for log 7V(Hi) = 15 
(thin lines) and logiV(Hl)= 17 (thick lines), [X/H] = -0.5 
(solid lines), [X/H] = -1 (dashed lines), and [X/H] = -1.5 
(dotted lines). The panels have different ionising spectra 
and silicon abudances: a,e HM ionising spectrum and so- 
lar abundance pattern; b,f SSP spectrum and solar abun- 
dance pattern; c,g HM ionising spectrum and [Si/H] = 
0.1; d,h SSP spectrum and [Si/H] = 0.1 The superposed 
points are the observed Mgn/Civ and Cii/Civ column 
density ratios for component 12 of the system at z ~ 1.87 



tios are shown as black points. Three models are consis- 
tent with the data: an HM ionising spectrum with so- 
lar abundance pattern and metallicity [X/H] ~ —1 for 
logJV(Hl)= 15 or [X/H] = -1.5 for logiV(Hl)= 17 (pan- 
els a,e); an HM ionising spectrum with overabundance 
of silicon, logiV(Hi)= 15 and metallicity [X/H] = —1.5 
(panels c, h); an SSP spectrum with overabundance of 
silicon and a metallicity [X/H] ~ —0.5 (panels d, g). The 
Si hi /Si iv vs. Sin/Siiv plot favour the latter model (see 
Fig. 0). 

In conclusion, we find that the two components at z ~ 
1.8675 and at z ~ 1.87008 can be explained by a model 
with a local ionising stellar source and a gas metallicity of 
~ 1/3 solar. 



not as broad as those observed e.g. by Kirkman & Tytler 
(19990 

The high uncertainty in the H I column density for 
most of the components and the blending of many in- 
teresting absorption lines (refer to Sect. 2.3) makes an 
abundance analysis for this system difficult. A fit to the 
column densities in component 2, at z ~ 1.8675, shows 
that the SSP spectrum model is best suited with log U ~ 
-2.36 ± 0.04, [X/H] ~ -0.4 ± 0.1 and an overabundance 
of silicon, [Si/H] ~ 0.1. 

In the case of component 12 at z ~ 1.87008 for which 
we do not have a reliable determination of H I column 
density, we run a grid of photoionisation models vary- 
ing the main parameters. On the basis of what is found 
from the previous modelisation, we consider metallici- 
ties between ~ 1/3 and ~ 1/30 solar, solar abundance 
pattern or an overabundance of silicon, [Si/H] = 0.1. 
Fig. ^ shows the predicted Siiv/Civ ratio as a function 
of Mgn/Civ and Cii/Civ ratios for the ionising param- 
eter varying between —4 < \ogU < —1. Observed ra- 



4. The Siiv/Civ ratio 

It is well known that the Si iv/C iv ratio in high-ionisation 
systems depends critically on the strength of the He n ion- 
isation edge (54 eV) break. Direct observations of the He H 



Lya absorption in QSO spectra (Jakobsen et al. 1994; 


Davidson et al. 1996|; 


Hogan et al. 1997; Reimers et al. 


1997 


; Zheng et al. 1998 


) show a marked decrease of the 



opacity for z < 3. Detailed investigation of the Siiv/Civ 



ratio have been carried out by Songaila & Cowic (1996); 
Bavaglio et al. (1997] ); |Giroux fc Shull (1997[ ) and |Songaila| 



(1998) among others. They all agree on the fact that the 
above ratio increases strongly between z = 2 and z = 4 
with a possible discontinuity around z — 3. In particular, 
(Songaila (1998 ) collected a sample of Civ systems with 
-/V(Civ)> 5 x 10 12 cm" 2 , and obtained median values of 
0.043±g gJ| for all the systems below z = 3, and 0.15±0.04 
above z — 3 ( errors are 1 a com puted using the median 
sign method). Boksenberg (1995 ) questioned this result, 
based on an analysis of the redshift evolution of the ion 
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0.1 - 



This work 

A Kohler et al. (1999) 




0.01 



Fig. 11. The values of the Si iv/C iv column density ratio 
for the single components of the metal systems in our sam- 
ple (see text) are reported as a function of redshift. The 
two long-dashed lines are the median values computed by 
Songaila (1998) from lines in her sample below and above 
z = 3 



ratios in individual Voigt components within complex sys- 
tems. 

The Siiv/Civ ratios for the individual components of 
the three considered systems are plotted as a function of 
redshift in Fig. |ll] together with the data by Kohler ct al 



(1999|]]who observed 4 systems at similar resolution. All 
systems have z < 3. Data are spread over more than an 
order of magnitude and most of them are above the value 
estimated by Songaila (1998 ) for systems with z < 3. It 
has to be noted that all but one of the systems in our 
sample have total Civ column density > 10 14 cm -2 , while 
only 10 % of Songaila's systems at z < 3 belongs to this 
range of column densities. 

As a comparison, we consider the HST data on ab- 



sorpti on lin es arising in the Milky Way halo gas (Savage 
et al 



|000| ). The observed Si IV and ClV transitions are 
weak, thus we transform the observed equivalent width in 
column density by using the curve of growth relationship 
for optically thin gas. The weighted mean of the 9 ob- 
tained estimates (corresponding to different lines of sight) 
is Siiv/Civ ~ 0.413 ± 0.009 and the values are spread 
between ~ 0.1 and ~ 0.6. 

This result supports the idea that high column den- 
sity systems, which most probably arise in halos of galax- 
ies, have higher values of the Siiv/Civ ratio compared to 
weak systems arising in regions more typical of the inter- 
galactic medium. The systems studied here are probably 



influenced by local ionising sources rather than by the dif- 
fuse UV background. 

5. Conclusion 

We have discussed the nature of the three prominent metal 
absorption systems observed in the UVES spectrum of the 
HDF-South quasar J2233-606. 

We analyse the velocity structure of the systems 
and their ionisation and abundance properties mainly 
by means of photoionisati on models per formed with the 
Cloudy software package ( Ferland 1997 ). The main con- 
clusions we draw are the following: 

1. - The signature of a probable expanding flow in an 
Hn region is observed in the main metal component of 
the Lyman limit system, at z ~ 1.9259, which shows 
a correlation between ionisation state and kinematics. 
Furthermore, the two main components in the H I absorp- 
tion profile are shifted by —43 km s _1 and +184 km s _1 
with respect to that component. The observed column 
density pattern is explained by a photoionisation model 
with a local stellar source and three different metallici- 
ties, ~ 1/500 and ~ 1/100 solar for the two Hi clouds and 
~ 1/8 solar for the metal component. The high metallicity 
in the latter cloud and the marked inhomogeneity suggest 
that the line of sight intersects a region of ongoing star 
formation and recent metal enrichment. 

2. - The system at z ~ 1.94 shows high metallicity, ~ 1/4 
solar, when modelled with photoionisation by an HM spec- 
trum. There are three components at higher ionisation 
where strong O vi absorption is detected. The metallic- 
ity of these components is lower or equal to the one in 
the low ionisation components depending if the same or 
a different ionising source (SSP spectrum or collisional) 
is considered. Alii and Aim are badly underestimated by 
the model, this could be due to a still low ac curacy in the 
recombination coefficients for this element (Petitjean et 
|al. 1994Q . 

3. - The system at z ~ 1.87 shows associated Ovi ab- 
sorptions. The same ionisation origin for these lines and 
the narrow transitions of Civ, Si iv, Si ill, etc, is excluded, 
thus a more complex scenario has to be taken into account. 
In order to test the different explanations which have been 
proposed it will be interesting to obtain spectra at high 
resolution and signal-to-noise ratio of the O VI complexes. 
We could investigate only two of the 17 components, both 
agree with [X/H] ~ 1/3 solar and a local stellar ionising 
source. 

4. - The analysis of the Si iv/C iv column density ratio for 
a sample of systems at z < 2 with total iV(Civ)> 10 1 4 
cm -2 , shows a ditribution skewed towards higher values 
than observed by Songaila (1998] ). A comparison with 
analogous results for the gas in the halo of our galaxy im- 
plies that higher total C iv column density systems present 
higher values of the ratio Siiv/Civ, likely because they 
arise in a galactic environment and are influenced by local 
ionising sources rather than by the UV background. 
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Comp. 


z 


b 

(km s" 1 ) 


Ion 


logiV 

(cm- 2 ) 


1. . . 


1.86733 c 


60 ±6 


Civ 


13.31 ±0.09 


2. . . 


1.867494" 


5.3 ±0.2 


Si iv 


12.97 ±0.01 




(f) 


5.2 ±0.4 


Sim 


12.65 ±0.03 




(f) 


(f) 


Sin 


11.7 ±0.1 




(f) 


(f) 


Cn 


12.4 ±0.1 






5.9 ±0.3 


Civ 


13.75 ± 0.02 


3... 


1.86756 6 


25 ±2 


Civ 


13.63 ±0.05 


4... 


1.867748 a 


4.6 ±0.4 


Si iv 


12.49 ± 0.02 






5± 1 


Sim 


11.93 ±0.05 




1.867758° 


6. ±0.3 


Civ 


13.82 ±0.02 


5... 


1.867980 6 


1.5 ±0.7 


Civ 


12. ±0.2 


6... 


1.868459 6 


3±3 


Si iv 


11.6 ± 0.1 






9 ± 2 


Civ 


12.56 ±0.05 


7... 


1.868749 ft 


7.8 ±0.9 


Civ 


12.77 ±0.04 


8... 


1.868966" 


8.4 ±0.5 


Civ 


13.12 ±0.02 


9... 


1. 869414" 


1.5 ±0.8 


Civ 


12.5 ±0.2 


10... 


1.86964 c 


23 ±4 


Si iv 


12.3 ±0.1 






18 ± 1 


Civ 


13.62 ±0.04 


11. . . 


1.869767" 


6.6 ±0.2 


Si iv 


13.12 ±0.02 




(f) 


8.4 ±0.5 


Sim 


13.09 ±0.05 




1.869756 6 


8± 1 


Mgn 


11.91 ±0.04 








Aln 


11.1 ± 0.1 








Aim 


< 11.48 ±0.09 




1.869779" 


7.7 ±0.5 


Civ 


13.55 ±0.04 


12. . . 


1.870075 a 


8.3 ±0.6 


Si iv 


12.72 ±0.02 




(f) 


6.5 ±0.5 


Sim 


13.01 ±0.07 




1.870086 a 


4.8 ±0.6 


Sin 


12.08 ±0.09 








Mgn 


11.97 ±0.03 








Aln 


11 ±0.1 








Aim 


< 11.2 ±0.1 






10 ±2 


Cn 


13.53 ±0.04 




1.870061" 


12 ± 1 


Civ 


13.33 ±0.07 


13... 


1.870232 6 


4± 1 


Si iv 


12.01 ±0.08 




(f) 


6± 1 


Sim 


12.13 ±0.06 






33 ±6 


Civ 


13.1 ±0.1 


14. . . 


1.870411 fc 


3±2 


Si iv 


11.6 ± 0.1 




(f) 


5±1 


Sim 


11.7 ±0.1 


15... 


1.870582" 


4.2 ±0.5 


Si iv 


12.24 ±0.03 




(f) 




Sim 


12.17 ±0.07 






6.1 ±0.7 


Civ 


12.94 ±0.05 


16... 


1.871079" 


8.6 ±0.2 


Si iv 


12.86 ±0.01 




(f) 




Sim 


12.64 ±0.07 




(f) 


6.6 ±0.8 


Cm 


15.9 ±0.5 




1.871102" 


10.2 ±0.2 


Civ 


13.68 ±0.01 


17... 


1.871796 6 


6± 1 


Si iv 


12.03 ±0.05 




(f) 




Sim 


11.6 ±0.1 




(f) 


11 ±3 


Cm 


13.1 ±0.1 




1.871808" 


8.5 ±0.3 


Civ 


13.28 ±0.01 


Error < 5 x 1(T 6 , b 


Error 6^15 x 10~ 6 , c 


Error > 50 x 10" 



